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CD 1 Abstract. The double-exchange model, which has been extensively studied in the context of 

5_l ■ colossal magneto-resistance in perovskite manganese oxides, is known to exhibit a ferromagnetic 

' metallic state at low temperatures because of the interplay between localized moments and 

j t , itinerant electrons through the Hund's-rule coupling. Here we investigate numerically the 

■ ferromagnetic transition in the double-exchange model defined on the frustrated pyrochlore 

lattice as a simple model for ferromagnetic pyrochlore oxides. We demonstrate that the 
I . finite-size corrections are largely reduced by implementing averages over the twisted boundary 

' O ' conditions in the Monte Carlo simulation, which enables to estimate the ferromagnetic transition 

temperature in relatively small size clusters. The estimate is compared with that for the non- 
frustrated cubic lattice system. 

> : 

^ ■ 1. Introduction 

\^ \ The double-exchange (DE) model is a minimal model which explicitly incorporates the interplay 

between itinerant electrons and localized magnetic moments. The Hamiltonian is given by 

(N ■ 

O; ~H = -^2 t(c} a c ja + h.c.) - Jn^2 Si -Si - rii, (1) 

(ij)a i i 

where C{ a (c\ a ) is an annihilation (creation) operator of an electron with spin a at site i, t is 
the transfer integral for nearest-neighbor sites (ij), Sj and Sj are spin operators of itinerant 
electron and localized spin, respectively, which are coupled by the Hund's-rule coupling Jh, 
n i = J2a c lcr c i& i s the density operator, and fi is the chemical potential. The localized spins 
Sj are treated as classical vectors. The model was originally introduced by Zener pQ, and 
has been studied for understanding the physics of perovskite manganese oxides [2, 3J. In 
particular, the rediscovery of colossal magneto-resistance phenomena have stimulated extensive 
studies for the DE systems, including some extensions of the model such as the super-exchange 
interaction between localized moments, orbital degeneracy of itinerant electrons, and electron- 
phonon couplings [H [5] . 

The model (JJ is known to exhibit a ferromagnetic metallic state at low temperature (T) 
to gain the kinetic energy of electrons by aligning localized moments in parallel. This is called 
the DE mechanism [TJ. There have been many efforts to estimate the ferromagnetic transition 
temperature T c [6]. For the model on the three-dimensional cubic lattice, T c was determined 



precisely by a large-scale Monte Carlo simulation: For example, T c /t = 0.136(2) at the electron 
density n = 0.5 in the limit of Jn/t — > oo [7]. 

In this contribution, we present our numerical results for the ferromagnetic transition when 
the model is defined on the geometrically frustrated lattice structure. Among many frustrated 
lattice structures, we consider the pyrochlore lattice, which is a three-dimensional network of 
corner-sharing tetrahedra as shown in Fig. [D^a). One of the experimental motivations is found 
in a family of Mo pyrochlore oxides -R2M02O7: When the ionic size of rare earth element R is 
relatively large such as R = Nd and Sm, the compounds become ferromagnetic metal at low T 
[51 E] , and it was pointed out by the first principle calculations that the DE mechanism plays 
a key role for this behavior [10] , As a first step toward the understanding of thermodynamic 
properties of the pyrochlore systems, below we will investigate the ferromagnetic transition by 
Monte Carlo calculations for the simplest case, i.e., the model ([T]) in the limit of Jti/t — > 00. 




Figure 1. (a) Pyrochlore lattice structure. The box indicates the cubic unit cell, (b) Density 
of states per site for the non-interacting model on the pyrochlore lattice. 



2. Monte Carlo Simulation with Averaging over Twisted Boundary Conditions 

We employ a Monte Carlo method to take account of large fluctuations in the frustrated system. 
The method is a standard one in which configurations of classical localized spins are sampled by 
Monte Carlo procedure; the Monte Carlo weight is calculated by the exact diagonalization of the 
fermion Hamiltonian matrix for a given spin configuration. The bottleneck of the calculations 
is the exact diagonalization, which usually limits the accessible system sizes to several hundreds 
sites. An improved method based on the polynomial expansion of the density of states was 
developed by the authors I12j. but for frustrated systems, in general, the method becomes 
less efficient because of singular form of the density of states: Indeed, in the present pyrochlore 
case, the density of states exhibits a 5- functional peak due to two flat bands (co/t = 2) as well as 
two van-Hove singularities {oo/t = and —4) in the non-interacting model as shown in Fig. [2(b) . 
Therefore we here employ the standard method. 

To enable systematic analysis within the limited system sizes, we apply a technique of 
averaging over the twisted boundary conditions |13^. 114] . In this technique, a twisted boundary 
condition is imposed with replacing the transfer integral t by texp(icp ■ Sij), where </> denotes 
a magnetic flux and Sij represents the vector connecting the nearest-neighbor sites (ij). The 
average is taken by the integral over <fi, which is approximately calculated by the sum over 
grid points. It has been shown that the procedure reduces finite size effects originating from the 
discreteness of the wave numbers. 

It is shown for the model ([T]) that in the limit of — > 00, the averaging procedure for 
a state of a finite size cluster with a given set of spin configurations provides results for the 
infinite size system which consists of a periodic array of the finite size cluster. For example, 
for the perfectly ordered ferromagnetic state, the averaging procedure with — > 00 for any 



finite size cluster gives the exact result in the thermodynamic limit. In order to take this 
advantage, we apply the averaging procedure to each Monte Carlo snapshot for calculating the 
Monte Carlo weight. This corresponds to ensemble average over the independent systems with 
different boundary conditions. In the following calculations, we take the average over the grid 
points with <f> = ((^ m x — 1)^^, (2m y — l)A<fi y , (2m z — l)A(p z ), where A<p v = ir/2L v l v and 
m u = 1, 2, ■ ■ ■ , l v (y = x,y, z). Here, N<p = l x x ly x lz and L v is a linear dimension of the system 
measured in the cubic unit cell, i.e., the total number of sites N s = L x x L y x L z x 16. 

We demonstrate here the efficiency of the averaging technique by calculating the electron 
density n as a function of T. We show the Monte Carlo results at [i = in Figs. [2] and as an 
example. In Fig. [21 we present dependence for two different system sizes. In both cases, the 
results converge onto a single curve as increasing N^. Furthermore, the necessary Na, for the 
convergence becomes smaller as the system size iV s increases. Figure shows N s dependence 
with and without taking averages. In the case without taking averages, the results are largely 
scattered for different system sizes as shown in Fig. [3](a) . In contrast, as in Fig. E^b), when we 
take averages over a sufficient number of Na, chosen for each system size, all the results converge 
onto a single curve and the finite size effects are sufficiently small, even for the rectangular- 
shaped clusters. We have checked the efficiency in wide parameter regions of /i and T and 
confirmed that the averaging procedure is efficient enough to suppress the finite size effects 
coming from the discreteness of the wave numbers. We adopt used in Fig. E^b) for each N s 
in the following calculations. 
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Figure 2. T dependence of the electron density n at fj, = for different numbers of grid 
points N^. (a) The system size iV s = lxlxlxl6 and (b) N g = 2 x 2 x 2 x 16. The lines are 
guides for the eyes. 




Figure 3. T dependence of n at \i = for different system sizes iV s (a) without taking 
average {N& = 1 x lxl) and (b) with taking average Na, (the numbers are shown in the 
legend). The lines are guides for the eyes. 



3. Results and Discussion 

Applying the method above, we investigate the magnetic behavior of the model (pQ) in the limit of 
J~a/t — > oo at the electron density n = 0.5. The results are presented in Fig. [4j At T/t ~ 0.14, 
the square of total magnetization per site m = J2iSi/N s grows rapidly [Fig. [3(a)] and the 
uniform magnetic susceptibility \ = (d m | 2 ) — ( | m | ) 2 ) ^V s /X 1 exhibits a peak [Fig. 0(b)], signaling 
a ferromagnetic transition. We estimate T c from the crossing point of the Binder parameter 
g = 1 — ((|m| 2 ) 2 )/3(|m| 2 ) 2 for different system sizes: The estimate is T c /t = 0.135(15). 
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Figure 4. T dependences of (a) square of total magnetization per site, (b) uniform magnetic 
susceptibility, and (c) Binder parameter for the total magnetization per site. The electron 
density is fixed at n = 0.5 by controlling \i. 

The value of T c is quite similar to that for the cubic lattice model, T c /t = 0.136(2) [7j. This 
is reasonable because the DE ferromagnetism is governed by the kinetic energy of electrons, and 
does not strongly depend on the details of the lattice structure. 

Our results demonstrate that systematic and quantitative study of the phase diagram is 
feasible within relatively small size clusters by implementing the averaging procedure over 
twisted boundary conditions in the Monte Carlo calculations. This gives a starting point for 
further study of frustrated DE-based models and for understanding the physics of pyrochlore- 
based ferromagnets such as Mo pyrochlore oxides. Studies of phase competitions in the model 
including the super-exchange coupling between localized moments are in progress. 
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